Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN UT^X style file vl.4) 



Effects of quintessence on observations of Type la 
SuperNovae in the clumpy Universe 

M. Sereno/'20 E. Piedipalumbo,^'^ M.V. Sazhin^'^ 

^ Dipartimento di Scienze Fisiche, Universita degli Studi di Napoli "Federico 11", 

Comp. Univ. Monte S. Angela, Ed. G, Via Cinthia, 80126, Napoli, Italia. 

^ Istituto Nazionale di Fisica Nucleare, sez. Napoli, 
... Comp. Univ. Monte S. Angelo, Ed. G, Via Cinthia, 80126, Napoli, Italia. 

' Sternberg Astronomical Institute, Moscow, Russia. 
C " * Osservatorio Astronomico di Capodimonte, Via Moiariello, 16, 80131 Napoli, Italia. 

o 

2 February 2008 

CD 

^ . 

ABSTRACT 

We discuss the amplification dispersion in the observed luminosity of standard candles, 
like supernovae (SNe) of type la, induced by gravitational lensing in a Universe with 
dark energy (quintessence) . We derive the main features of the magnification proba- 
bility distribution function (pdf ) of SNe in the framework of on average Friedmann- 
' Lemaitre- Robertson- Walker (FLRW) models for both lensing by large-scale structures 

and compact objects. Analytic expressions, in terms of hypergeometric functions, for 
On ' luminosity distance-redshift relations in a flat Universe with homogeneous dark energy 

. have been corrected for the effects of inhomogeneities in the pressureless dark matter 

I (DM) . The magnification pdf is strongly dependent on the equation of state, wq , of the 

. quintessence. With no regard to the nature of DM (microscopic or macroscopic), the 

' dispersion increases with the redshift of the source and is maximum for dark energy 

with very large negative pressure; the effects of gravitational lensing on the magnifi- 



O ' cation pdf, i.e. the mode biased towards de-amplified values and the long tail towards 

^ , large magnifications, are reduced for both microscopic DM and quintessence with an 

^ ' intermediate wq. Different equations of state of the dark energy can deeply change 

the dispersion in amplification for the projected observed samples of SNe la by future 
space-born missions. The "noise" in the Hubble diagram due to gravitational lensing 
, strongly affects the determination of the cosmological parameters from SNe data. The 

$H ' errors on the pressureless matter density parameter, ^Im, and on wq are maximum for 

quintessence with not very negative pressure. The effect of the gravitational lensing 
is of the same order of the other systematics affecting observations of SNe la. Due to 
the lensing by large-scale structures, in a flat Universe with JIm = 0.4, at 2 = 1 a 
cosmological constant {wq = —1) can be interpreted as dark energy with wq < —0.84 
(at 2-(T confidence limit). 

Key viTords: cosmology: theory - dark matter - distance scale - gravitational lensing 
- large-scale structure of the Universe - supernovae: general 



1 INTRODUCTION 



During the last years, two independent groups, the High-z SuperNova Search Team (Schmidt et al. 1998) and the Supernova 
Cosmology Project (|Perlmutter et al. 199S: ) have given strong evidences of the acceleration of the Universe's expansion ( Riess 



et al. 1998; 



Perlmutter et al. 1999). Several other observational and theoretical evidences, like measurements of the anisotropy 



of the Cosmic Microwave Background Radiation ( ie Bernardis et al. 2000 ) and inflationary theories, strongly support a 
flat or nearly flat Universe. On the other hand, direct measurements of SIa/ from dynamical estimates or X-ray and lensing 
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observations of clusters of galaxies indicates that Qm is significantly less than unity, Qm — 0.3 ( Turner 2000| ). To solve this 
puzzle, a new type of energy component in the Universe, now called dark energy or quintessence, was proposed; dark energy 
with strongly negative pressure is required to explain acceleration {wq = pg/pq < —1/3, where pg and pq are, respectively, 
the pressure and energy density of the dark energy). 

Observations of SNe la, which at low redshifts are sensitive to the deceleration parameter qo = {^Im + (1 + 3wQ)nQ)/2 
where fig is the density parameter of the dark energy, rely on several properties of these sources. SNe la are very luminous 
and have a small intrinsic dispersion in their peak absolute magnitude, SM ~ 0.3 (Filippenko & Riess 200C). These features 
make them the long expected standard candles for cosmology. 

A standard candle is a source with known intrinsic luminosity (L) (or absolute magnitude). Measurements of its apparent 
flux {J^) allow us to determine the photometric distance Dl to the source via equation 



Dr 



L 



(1) 



Using standard candles, it is possible to plot the Hubble diagram (that is, the redshift of an object versus cosmological distance 
to it or vice versa) with very high precision and estimate the global cosmological parameters. 

There are several candidates for the dark energy. The oldest one, initially introduced by Albert Einstein as a new funda- 
mental constant of nature, is the cosmological constant (wg = — 1). After the formulation of inflationary theory, cosmologists 



found that a A term can be introduced dynamically (Dolgov, Sazhin & Zeldovich 1990; Zel'dovich 1992; Sahni & Starobinski 



200C 
198^ 


) 


a dynamical A term by a scalar field slowly rolling down its potential (wg > — 1) ( 


Peebles & Ratra 1988 




Wetternich 




□striker & Steinhardt 1995 


; Caldwell et al. 1998; Zlatev, Wang & Steinhardt 1998; de Ritis et al. 2000; 


Rubakov 200C ; 


Rubano & Scudellaro 2001 


) can support a static energy component with positive energy density but negative pressure. Other 



poss ibilities for the quintes sence are represented either by networks of light, non-intercommuting topological defects ( Vilenkin 
1984 ; ^pergel fc Pen 1997 ) (mq = — m/3, where m is the dimension of the defect: for a string wg — — 1/3; for a domain 
wall Wg = —2/3) or by the so called X-matter (Chiba, Sugiyama & Nakamura 1997; Turner & White 1997). Alternatively to 



quintessence, a Universe filled with Chaplygin gas (Kamenshchik, Moschella & Pasquier 2001) is an additional alternative to 
obtain a negative pressure. Generally, the equation of state wg evolves with the redshift, and the feasibility of reconstructing 
its time evolution has b een investigated ([Cooray fc Ifuterer 1999 ; Chiba & Nakamura 200C 
200l |; I Sahni et al 20o6|; [Goliath et al. 200l|; [Huterer fc Turner 2 00l|; [Nakamura fc Chiba 2001 



Maor, Brustein & Steinhardt 



Pavlov et al. 2002; Wang fc 



Garm vich 2001; Wang fc Lovelace 2001; WeUer fc Albrecht 2001s ; Weller fc Albrecht 2001b Vamamoto fc Futamase 2001 



Gerke & Efstathiou 2002). Since in flat FLRW models the distance depends on wg only through a triple integral on the 



redshift (Maor, Brustein & Steinhardt 2001), wg(^z) can be determined only provided a prior knowledge of the matter density 
of the Universe (Goliath et al. 2001). In what follows, we will consider only the case of a constant equation of state. 

Astrophysical sources other than SNe have been long investigated to build the Hubble diagram. Two independent lumi- 



Reichart et al 



nosity estimators, the first one based on the variability of Gamma-Ray Bursts (GRBs) (Reichart & Lamb 2001 
200l|)]and the second one derived from the time lag between peaks in hard and soft energies ([Norris, Marani fc Bonnell 2000[), 



have been recently proposed to infer the luminosity distance to these sources. On the other hand, standard rods, as compact 



radio sources (Gurvits et al. 1999) or double radio galaxies (GJurevich & Zybin 2000), have been long studied to evaluate the 



angular diameter distance to cosmological sources. With no regard to their different physical origins, all these observations 
are affected by gravitational lensing of the sources. In this paper, we want to quantify the effect of inhomogeneities in the 
pressureless matter on the determination of the distance. In particular, we will study SNe la, whose importance in the deter- 
mination of the cosmological parameters makes necessary a complete study of all systematics. For light beams propagating in 
the inhomogeneous Universe, the expression of the luminosity distance in terms of the cosmological parameters, as obtained 
from Eq. (fH), changes with respect to the corresponding FLRW model. Mathematical considerations for non-flat models of 



Universe are done in Sereno et al. (2001). In this paper, we will discuss the simple case of the flat Universe and the influence 
of dumpiness on the Hubble diagram. We consider inho mogen eous pressureless matter and smooth dark energy. For the more 
general case of inhomogeneous quintessence see Linder (1988) and Sereno et al. (2001). 

The paper is as follows. In Sect. 2, we introduce the on average FLRW models; we discuss the Dyer-Roeder (DR) equation 
and its analytical solution in terms of hypergeometric functions. In Sect. 3, we consider the case of the homogeneous Universe. 
Section 4 contains the discussion on the statistical nature of the lensing dispersion induced either by large-scale structures or 
compact objects; we study the magnification pdf induced by gravitational lensing and the connected systematic errors in the 
estimate of and wg. Section 5 is devoted to some final considerations. 



2 THE LUMINOSITY DISTANCE-REDSHIFT RELATION 

The standard Hubble diagram is computed with relations that hold in FLRW models, that is assuming all gravitating 
pressureless matter homogeneously distributed. Instead, observational data are taken in the inhomogeneous Universe, where 
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sources most likely appears to be de-magnified relative to the standard Hubble diagram. In fact, light bundles, propagating 
far from clumps along the line of sight from the source to the observer, experience a matter density less than the average 
matter density of the Universe. Several attempts have addressed the problem of the redshift dependence of the distance in a 
clumpy Universe: by relaxing the hypothesis of homogeneity and using the Tolman-Bondi metric instead of the FLRW one 



(Celerier 200C); by quantifying the small deviations from the isotropy and homogeneity of the Ricci scalar (Trentham 2001) 



by considering a local void (Tomita 2001). One of the historically more important and widely used framework is the on-average 



FLRW Universe (Schneider et al. 1992; Seitz, Schneider & Ehlers 1994), where it is assumed that the relations on a large scale 



are the same of the corresponding FLRW Universe, while inhomogeneities only affect local phenomena like the propagation 
of light. A simple way to account for this scenario is to introduce the smoothness parameter a a/, that, in a phenomenological 
way, represents the magnification effect experienced by the light beam (Tomita 1998; Wang 199t). In general, the smoothness 
parameter is redshift dependent since the degree of smoothness in pressureless matter evolves with time. When oa/ = 1, 
we reduce to the FLRW case; when aM < 1, there is a defocusing effect; om ~ represents a totally clumped Universe. 
This limiting case, sometimes known as "empty-beam approximation", corresponds to the maximum distance to a source 
for light bundles which have not passed through a caustic (Schneider et al. 1992). Historically, qa/ is defined as the fraction 
of pressureless matter smoothly distributed. The distance recovered in on average FLRW models, sometimes known as DR 



distance, has been long studied (Zel'dovich 1964 
198d 



Linder 1988 



Kantowski 196S: 



Dyerfc Roeder 1972; Dyer fc Roeder 197; ; Korolev fc Sazhin 



; Beitz, Schneider & Ehlers 1994) and now is becoming established as a very useful tool for the interpretation 
of experimental data (Kantowski 1998, Kantowski & Thomas 2001; Perlmutter et al. 1999; Giovi & Amendola 2001). The DR 



equation for the angular diameter distance. Da, in the case of a Universe with inhomogeneous pressureless matter and dark 
ergy with constant equation of state, has already been obtained starting from the optical scalar equations ( Linder 1985 ; 
ovi fc Amendola 200 1[ ) or using the multiple lens-plane theory (Sereno et al. 2001). For flat Universes, Qm + — 1, it is 

(1 + zf (1 + Ml + ^r'') ^ + ^ (7 + {3wQ + 7M1 + .)«-«) 



dDA 



dz 



(2) 



+ -{aM + {wQ + l)ii{l + zf 



•Da = 0, 



where n = 75-^ and aM is the smoothness parameter. The boundary conditions on Eq. (|2[) are 



dz 



Da{0) 



Da 



0, 
c 



(3) 



It is straightforward to obtain the corresponding equation for the luminosity distance. Using the Etherington principle ( ^ilther- 



ingtoi 1933), 



Dl = {l + zfDA, 

we can substitute in Eq. (Q), 



il + zY{l+ti{l + z) 



+ 



(^ + ^)(l + (l-3u,g)Kl + .)^™«)'''" 



' dz'^ 2 

flil + zf^Q^DL 



dz 



(4) 



(5) 



V 2 2 
the boundary conditions are, again. 



Dl{0) = 0, 



(6) 



—D 

dz ^ 



c 



The solution of Eq. satisfying the boundary conditions in Eq. (6), takes the form 
c Di(0)D2(z)-Di(z)-D2(0) 



Dl[z) = 



Ha 



W{Q) 



where Di{z) and 1)2(2) are two linearly independent solutions of Eq. (p|) and W{z) = Di{z 



D2M dDi(z) 



dz 



dz 



(7) 

D2{z) is the 

Wronskian of the solutions system. In what follows, we will consider the DR equation with a constant smoothness parameter. 
We introduce the parameter 

V25 - 24QAf 



/3 



To solve Eq. (^), we perform the transformation of both the independent and dependent variables, 
u = -fi{l + zf'''3, Dl{z) = u'''-Q Rl{z). 



(8) 



(9) 
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Figure 1. The luminosity distance in a flat and smooth Universe with f2]\/ = 0.3. The distance is in units of c/Hq. 



With such a transformation, Eq. Q reduces to the hypergeometric equation for Rl, 



cPRl 



+ 



1 + 



cLRl 
du 



3wQ I u 2(1 — It) 
A pair of independent solutions of Eq. (|1C|) is 

"4/3-14/3 + 1 1 2/3 



12wQ 



4/3 -1\/^ 4/3 + 1^1 



l2wQ 2 / u{l-u) 



-Rl = 0. 



Ri{u) 



2F1 



i?2(u) ^ 2F1 



12wQ ' 12wQ 2 ' 3wQ 

4/3 + 1 -4/3 + 11 2/3 



(10) 



(11) 



+ -, hi, It 

12wQ 12wQ 2 3wQ 



where 2^1 is the hypergeometric function of the second type. Inserting the expressions for Ri and R2 in Eq. (g) and substituting 
in Eq. (Q), we have the final expression for the luminosity distance, 

c 1 



Dl{z) 



Ho 2/3^n^ 

X <j (1 + 2)^+" 2F1 



(12) 



4/3 + 1 1 1-4/3 ^ 
12wQ ' 2 12wQ ' 



2/3 



X2F1 



4/3-1 1 , 4/3 + 1^^ ^ ^2/3 ^ ^M-l ^^ , ^,^3^, 



12wQ 2 12wQ 



3iOQ SIm 



2F1 



4/3+1 1 ^ 1-4/3 ^ 2/3 - 1 ^ .3 

12toq ' 2 12mQ ' 3iOQ ' SIm 



X 2F1 



4/3-1 1 ^ 4/3 + 1 2/3 ^M-l 



12WQ 2 12WQ 



3wQ f2j\ 



For the case of a cosmological constant, wq = —1, Eq. (12) reduces to equation (16) in Kantowski & Thomas (2001), as 
we can see by using the property of the hypergeometric functions 

1 



2F1 [a, b, c, x] 



(1 - xY 



2F1 



a,c ~ b, c, 



x-1 



(13) 



and noting that the clumping parameter v in Kantowski & Thomas (2001) corresponds to (/3— 1)/2. The case of the cosmological 
constant is also studied in (Kantowski 1996; Kantowski, Kao & Thomas 200C; Dcmianski et al. 2000| ). 



3 THE LUMINOSITY DISTANCE IN THE HOMOGENEOUS UNIVERSE 

For the FLRW case {aM = 1), Eq. (B is solved by 



Dl{z) = -^{1 + z) ^ 

^0 Jo ^nM(i + z'y + (1 - nM)(i + 2')«(»«+i) 

This expression is equivalent to Eq. (12) when /3 — 1/4, 

nM - 1 



dz 



Ho vJTj 



- — ,i 1-— . 

6wQ 2 6wQ 



(14) 



(15) 
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2^1 



&WQ 2 



1 ,!lH_i(l + ,)3», 



Our Eq. ( p^ is equivalent to the expression found by Bloomfield Torres & Waga ( 1996| )(see also Giovi & Amendola (2001)). 
In the Einstein-de Sitter case {Qm = 1 or wq 0), Eq. ( p^ reduces to 

(16) 



i3i(z) = — 2(l + z)( 1- 



as can also be seen directly by solving the integral in Eq. ( p^ ) . Figure (0) plots the luminosity distance for different equations 
of state: the distance increases for decreasing wq. 



4 THE MAGNIFICATION PROBABILITY DISTRIBUTION FUNCTION 



The amplification of a source at a given redshift has a statistical nature. For narrow light-beams, the effect of gravitational 
lensing results in the appearance of shear and convergence in images of distant sources according to the different amount and 
distribution of matter along different lines of sight. So, gravitational lensing increases the level of errors in the Hubble diagram 
([Kantowski et al. 1995|; [Kantowski 1994 ^rioman 1997[ [Wambsganss et al. 1997| ; [HoIz 1998| ; [HoIz fc Wald 1998| [Metcaff 199E| ; 



Barber 200C; Porciani & Madau 2000). In the framework of the on average FLRW Universe, we can account for this effect by 



considering a direction dependent smoothness parameter om. Now, om represents the effective fraction of matter density in 



the beam connecting the observer and the source and depends on the distribution of matter in the beam (Wang 1999); values 
of om greater than 1 account for amplification effects. 

There is an unique mapping between the magnification ^ of a standard candle at redshift z and the direction-dependent 
smoothness parameter at z (Wang 1999). According to Eq. (|l|), the magnification fj, of the source with respect to the maximum 
empty-beam case (aj\/ = 0) is 



DLjctM = 0) 
DL{aM) 



(17) 



Once derived the magnification (that is, once found the distance by integrating the null-geodesic equation or using ray-tracing 
techniques along a line of sight) of a source at epoch z, the corresponding smoothness parameter is determined in comparison 
with the DR distance: the solution of the DR equation for that constant value of aM matches, at redshift z, that given value 
of the distance ( [Tomita 199*^ ; |Wang 1999| ). 

The shape of the magnification pdf depends on the redshift of the source, on the cosmological parameters and on the 



nature of the DM. The dark matter can be classified according to its clustering properties (Metcalf & Silk 1999; Seljak & Holz 



199£ 



Mortsell et al. 2001): microscopic DM consists of weakly interacting massive particles (WIMPs), such as neutralinos 



(Gurevich & Zybin 1995; Surevich, Zybin & Sirota 1997) and clumps on galaxy halo-scales; macroscopic DM consists of 
compact objects, such as massive compact halo objects (MACHOs) or primordial black holes. 

According to N-body simulations of large scale structures in cold dark matter models, galactic halos are expected to 
contain a large number of small substructures besides their overall profile. However, this type of small-scale structure does 



not act as a compact object and only clumps of galaxy-size contribute appreciably to the lensing (Mortsell et al. 2001) 



In the framework of the on average FLRW models, the fi-pdi is characterized by some general features with no regard to 
the nature of the DM. Under the assumption that the area of a sphere at redshift z centred on the observer is not affected by 
the mass distribution, the photon number conservation implies that the mean apparent magnitude of a source at z is identical 



to the FLRW value (Weinberg 1976; Schneider et al. 1992), 



{fl) = flFL = fJ-{aM = 1) > 1. 



(18) 



Since the matter is clumped, most of the narrow light-beams from distant sources do not intersect any matter along the line 
of sight resulting in a dimming of the image with respect to the fiUed-beam case: the mode of the pdf, fipeak , is biased towards 
the empty beam value. 



f^peak < (/i). 



(19) 



The third feature is a tail towards large amplifications which preserves the mean. In terms of the magnification relative to 
the mean. 



5^ ; 



DLjaM = 0) 



DL{aM = 0) 



DL{aM = 1) 



(20) 



the 5/i-pdf has the mean at Sfi — 0, the peak value at (5/ipeafe < and a long tail for 5fi > 0, with no regard to the source 
redshift and to the cosmological parameters. It follows from these very general considerations that a simple way to characterize 



© 0000 RAS, MNRAS 000, 000-000 



6 M. Sereno, E. Piedipalumbo, M. V. Sazhin 



the pdf is to consider the parameter A/i, defined as the difference in ampiification between the mean FLRW vaiue and the 
magnification in the empty beam case (qj/ — 0), A/i = —5^{aM ~ 0). When A^ increases, the mode vaiue moves towards 
greater de-magnification: to preserve the total probability and the mean value, the pdf must both reduce its maximum and 



enlarge its high amplification tail. From the properties of the angular diameter distance in a clumpy Universe ( sereno et al, 
200lj)] it follows that A^ increases with the redshift of the source and with dark energy with large negative pressure. So, 
the dispersion in the /i-pdf due to gravitational lensing increases with z and it is maximum for the case of the cosmological 
constant (see also Bergstrom et al. (2000)): quintessence with wq > — 1 reduces the bias towards large de-amplifications of 
the peak value of the pdf, partially attenuating the effect of the dumpiness. 



4.1 Lensing by microscopic dark matter 



The gravitational lensing effect by large-scale structure on the apparent luminosity of distant sources in the Universe has 



been studied either with N-body simulations ( Wambsganss et al. 1997 



Tomita 1998 



199e) 



Barber 2QQC: Barber et al. 2000: Jain et 



al. 2000) or with the integration of the geodesic deviation equation (Holz & Wald 1998; Bergstrom et al. 200C) in a Universe 
filled with either isothermal spheres or Navarro- Frenk- White profiles ( Navarro, Frenk fc White 1995| ; [Slavarro, Frenk fc White 



The /i-pdf for the smoothly distributed DM is characterized by two main trends with increasing redshift: an increase in 
the dispersion and an increasing gaussianity. As we look back to earlier times, the Universe becomes smoother on average 
and lines of sight become more filled in with matter: light bundles intersect more independent regions along their paths and 
the resulting /i-pdf approaches a gaussian by the central limit theorem ( ^eljak fc Holz 1999 ; Wang 1999). The corresponding 



QM-pdf also becomes symmetric but it reduces its dispersion and its mode goes to the fiUed-beam value ( Tomita 1998 ; Wang 



199£ 



Barber et al. 2000). The trends in dispersions in the /t-pdf and ctM-pdf are opposite since, with increasing redshift, a 



large variation in the distance corresponds to a small variation in the smoothness parameter (Sereno et al. 2001). Wambsganss 



et al. (1997) used the ray-tracing method for large-scale simulations in a cold dark matter Universe, normalized to the first 



year COBE data with Qm ~ 0.4, Q,q — 0.6, wq — — 1, with a spatial resolution on small scales of the order of the size of a 
halo, to derive the /i-pdf at different redshifts. Wang ( 199!: ) was able to find empirical formulae for the fitting of the /i-pdf 
and of the corresponding QAf-pdf, 



, {cum, 



Da(qm =0) 

2/t3/2 



As noted by Tomita ( [l998| ), the angular diameter distance depends on qm linearly for < 
can approximate 

ODa 



doiM 



Da^cxm = 1)- DA{aM = Q). 



(21) 



' 5 and, with high precision, we 



(22) 



In Fig. (g), we plot the (5/i corresponding to the mode of the aAf-pdf (as plotted in figure (2b) in Wang (199!;)) as a function 
of the redshift: while the mode value of the aAf-pdf goes to the fiUed-beam value for increasing redshift, the variation in 
magnification with respect to the FLRW mean increases; that is, the bias increases with z. 

To study the role of the quintessence in the magnification dispersion of standard candles, we consider the same matter 



content, that is the same aAf-pdf (Wang 199!:), for different equations of state. Models with different cosmological parameters 
produce, in general, different QfA/-pdf, predictable by numerical simulations; but, to consider the influence of the dark energy 
on the /i-pdf, it suffices to use the same matter distribution in Eq. (pl|). This is equivalent to assume that the dependence on 
quintessence enters Eq. (^l|) through the angular diameter distances and that the effect on PaM is of the second order. So, 
for analytical convenience, we can use the same Pcm derived in Wang (1999) for several cosmological models with the same 
Q.M but different equations of state. In Fig. (^, the /i-pdf is plotted for two source redshifts and for two different equations 
of state: the /t-pdf becomes more and more symmetric with z and the dark energy reduces both the dispersion and the bias. 

The effect of gravitational lensing by large-scale structure affects significantly the determination of the cosmological 
parameters from observations of standard candles. Observed SNe la represent individual sources at each redshift and do not 
sample evenly the probability distribution: at a fixed redshift, we will observe the mode value of the distribution and not the 
mean one (Wambsganss et al. 1997; Barber 2000). For Q.m — 0.4, wq — —1 and z = 1, the mode is fipeak = 1-14 and the 
magnification values above and below which 97.5% of all of the lines of sight fall are fiiow = 1.11 and Hhigh = 1.28. This 
dispersion induces uncertainties in determining f^Ai and the equation of state. Assuming a flat Universe with cosmological 
constant, a Universe with Qm = 0.4 will be interpreted as a model with ~ 0.42l[J jj only because of the gravitational 
lensing noise. Here and in what follows, the error bars represent 2-a limits. With the constraint of SIai ~ 0.4, a cosmological 
constant might be interpreted as dark energy with wq < —0.84. 

For a flat Universe with SIa/ = 0.4 and wq = —0.5, at z = 1 it is /ipeafc ~ 1.11, niom ~ 1.09 and fihigh ~ 1.23. With the 
constraint wq — —0.5, we should estimate Qm ~ OAS'^g'^l^; assuming Qm = 0.4, it is wq = — 0.461!^q'24. 
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Figure 2. The magnification relative to the mean calculated for the peak value of the aj\,/-pdf as found in Wang (1999). It is Qf,i = 0.4, 
f2o = 0.6 and wn = —1. 
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Figure 3. The amplification pdf as a function of Sfi, the magnification relative to the mean, for microscopic DM. The sharply peaked 
line are for z = 0.5, the smoother ones for z = 2. Solid and dashed lines correspond, respectively, to wq = —1 and —1/2. It is Qm = 0.4 
and Qq = 0.6. Solid and dashed lines have the same matter distribution but different cosmological backgrounds. 



Although the lensing dispersion is reduced in a quintessence cosmology, the errors induced on the cosmological parameters 



increase. The reason is that in this models the luminosity distance is less sensitive to the cosmology (Sereno et al. 2001) 



4.2 Lensing by compact objects 



The eflect of gravitational lensing is maximum when the matter in the Universe consists of point masses (Holz & Wald 



1996); as seen above, this case is not included in the small-scale structures in the microscopic DM (Mortsell et al. 2001). The 



universal fraction of macroscopic DM is still unknown. Direct searches for MACHOs in the Milky Way have been performed by 
the MACHO and EROS collaborations through microlensing surveys. According to the MACHO group (Alcock et al. 2000), 
the most likely halo fraction in form of compact objects with a mass in the range 0.1 — 1 Mq is of about 20%; the EROS 
collaboration (Lasserre et al. 2000) has set a 95% confidence limit that objects less than 1 Mq contribute less than 40% of 
the dark halo. However, the average cosmological fraction in macroscopic DM could be significantly diflterent from these local 
estimates. 

The properties of the ^-pdf are essentially independent of both the mass spectrum of the lenses (this statement is 



strictly true for point sources (Schneider & Weiss 1988)) and the clustering properties of the point masses, provided that the 



clustering is spherically symmetric (Holz & Wald 1998). The dispersion in luminosity of standard candles is non-gaussian, 
sharply pea k ed at the empty be am value and has a long tail towards large magnifications falling as fi~'^ (Paczyhski 1986; 

' Wald 19*9^ , caused by small impact parameter lines of sight near the compact objects; so, its second 



Ranch 1991; Holz 



moment is logarithmically divergent and the law of large numbers fails: if strongly lensed ev ents are removed from the data 
sample, a bias will be introduced towards smaller apparent luminousities (Holz fc Wald 199S). 



A comparative analysis of the /x-pdf in the case of either microscopic DM or compact objects has put in evidence two 
main differences: the high magnification tail is larger for macroscopic DM and the mode of the distribution is nearer the 



average value in the case of lensing by large-scale structures ( jeljak fc Holz 199S: ; Mortsell et al. 2001 ) 



The /i-pdf in a Universe filled with a uniform comoving density of compact objects depends on a single parameter, the 
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Figure 4. The ojv/-pdf for macroscopic dark matter. Solid and dashed hnes correspond respectively to 2; = 0.8 and 1.5. It is Qf,i = 0.3, 
Qq = 0.7 and wq = —1. 



mean mae nification (/i) ( [Rauch 199l| ; ^eljak fc Holz 199E| ). Based on Monte-Carlo simulations, Rauch ([L991|) gives the fitting 
formula 



p(m) 



„fc(M-l) 



3/2 



(23) 



where the parameter b is related to the mean magnification by 6 = 247exp [—22.3(1 — (fi) ■ The approximation holds for 

(n)^^^'^ ~ 0.8, a condition verified up to 2: ~ 2 in a Universe with low matter density, with no regard to the equation of state 

WQ. 

The OM-pdf corresponding to the distribution in Eq. ( p^ ) is highly non-gaussian, see Fig. (^. The pdf decreases mono- 
tonically from the empty beam value to high values of the smoothness parameter. With increasing redshift, the QM-pdf tends 
to flatten and the probability for the fiUed-beam case and for high values of qm grows. 

In Fig. (^, we show the (5/i-pdf for two source redshifts and for two values of wq: quintessence with wq > — 1 reduces 
the effect of dumpiness. For z = 0.5, the variation in the distance modulus from the empty-beam case to the fiUed-beam one 
is 0.033(0.039) mag for wq = -l/2(-l); for 2 = 1, it is 0.108(0.138) mag for wq = -l/2(-l); for z = 1.5, it is 0.205(0.268) 
mag for wq — — 1/2(— 1). For 2 ~ 1, the bias towards the empty-beam value can be compared with the dispersion of 0.17 mag 
in the peak magnitudes of SNe la after the application of methods as the "multi-colour light curve" method (Riess, Press & 
Kirsh jier 1996| ). The effect of gravitational lensing is of the same order of magnitude as the other systematic uncertainties 
that limit the conclusions on the cosmological parameters based on SNe la Hubble diagram dFihppenko fc Riess 2000| ). The 
correlation between host galaxy type and both luminosity and light-curve shape of the source; interstellar extinction occurring 
in the host galaxy and the Milky Way; selection effects in the comparison of nearby and distant SNe; sample contamination 
by SNe that are not SNe la can produce changes as large as 0.1 mag in the measured luminosities of SNe la. 

The effect on the estimate of the cosmological parameters of gravitational lensing by a totally clumped model with only 
macroscopic DM is quite dramatic. For a source redshift of z = 1, a Universe with Qm ~ 0.3 and a cosmological constant can 
be interpreted as a model with Qm = 0.42 and wq — —1 or as one with Qm ~ 0.3 and wq = —0.71. These systematic errors 
increase in a quintessence cosmology with wq > —1. For z = 1, a Universe with Qm = 0.3 and wq = —2/3 will be interpreted 
as a model with Qm = 0.45 and wq = —2/3 or one with Qm ~ 0.3 and wq = —0.46. 



5 CONCLUSIONS 



Observations of SNe la are strongly affected by inhomogeneities in the Universe. For redshifts 2 ~ 1, the variation in the 
distance modulus from a standard flat FLRW model to a clumpy Universe with the same content of pressureless matter 
can be considerably greater than other systematic effects. The effect of amplification dispersion by gravitational lensing 
must be accurately considered. The prospects of future space-born missions, like the SuperNova Acceleration Probe (SNAP 



- Http://snap. lbl.gov) and the Next Generation Space Telescope, of determining properties of the dark energy have been 
discussed ( |Gohath et al. 200l| ; [Weller fc Albrecht 2001a| ; [Weller fc Albrecht 2001b| ; [Gerke fc Efstathiou 200^ ). According to 



these studies, SNAP data should only distinguish between a cosmological constant and quintessence with wq relatively far 
from —1. When SNe observations are combined with an independent estimate of Qm, for example from galaxy clustering 



( Verde et al. 2002 ), the degeneracies among the q uintessence models can be significantly reduced and some constraints on the 



time evolution of the equation of state can be put (Weller fc Albrecht 2001b; Gerke fc Efstathiou 2002 ) . However, these studies 



only consider measurement errors and intrinsic dispersion of the sources, neglecting the systematic and redshift dependent 
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Figure 5. The magnification pdf for macroscopic dark matter as a function of 5^, the magnification relative to the mean. SoHd and 
dashed fines correspond, respectively, to wq = —1 and —2/3. It is Q,m = 0.3, Oq = 0.7. Left panel: the source redshift is z = 1; right 
panel: it is z = 1.5 
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Figure 6. Amplification dispersion relative to the mean due to gravitational lensing by macroscopic DM for the projected 1-year SNAP 
sample. Thick and thin lines correspond, respectively, to wq = —1 and —1/2. It is ^Im = 0.3, Oq = 0.7. Intrinsic dispersion of SN 
luminosities is not considered. 



error induced by gravitational lensing. We have shown how, also assuming an exact knowledge of Q,m, in the redshift range 
covered by future missions a cosmological constant can be interpreted as dark energy with wq > —1. For SIm = 0.4 and 
z; = 1, a A constant may be interpreted as quintessence with wq < —0.84, only due to the lensing by large-scale structure. 
A fraction of DM in form of compact objects will make the situation even more dramatic. So, also with a prior knowledge of 
the remaining cosmological parameters, gravitational lensing can make the statements on the properties of dark energy based 
on SNe data significantly less certain. 

The effect of inhomogeneities dominates at high redshifts and should be one of the main systematics in attempting to 



build the Hubble diagram with GRBs (Norris, Marani & Bonnell 2000; Reichart . 
Deng 



Lamb 2001; Reichart et al. 2001; Schaefer 



Band 2001). The physical origin of GRBs is still uncertain, but recent observations suggest that they are related to 



the violent death of massive stars. Under the hypothesis that GRBs trace the global star formation history of the Universe, 



their assumed rate is strongly dependent on the expected evolution of the star formation rate with the redshift (Porciani 



Maiau 2001). While some scenarios prefer a redshift distribution of the GRB comoving rate peaked between z 

> 



1 and 



2 and out to very high redshift (Porciani & 



2, according to other ones the comoving rate remains roughly constant at z 

Madai 2001). Furthermore, the lack of strong lensing events in the fourth BATSE GRBs catalog (Holz, Miller 



Quashnock 

1998|y suggests that, at the 95% confidence level, the upper limit to the average redshift of GRBs is ~ 3 in a flat, low-matter 
density Universe with cosmological constant. According to these considerations, the effect of gravitational lensing would be 
really dominant in the Hubble diagram built with GRBs. 

As an example, we consider the GRB redshift distribution derived from a combined analysis of two independent luminosity 



indicators (Schaefer, Deng & Band 2001). Examining a sample of 112 GRBs from the BATSE catalog, Schaefer et al. (2001) 



found redshifts varying between 0.25 and 5.9 with a median of 1.5. At z = 1.5, gravitational lensing by large-scale structures, in 
a model with 57a/ = 0.4 and wq = —2/3, induces a magnification distribution with fipeak = 1.25, fiiom = 1.20 and fihigh ~ 1.46. 
Assuming wq = —2/3, we will estimate Qm = 0.43lQ'5g; assuming Q,m ~ 0.4, we will estimate wq < —0.51. 

Although the lensing dispersion on the luminosities of standard candles represents a noise in the determination of the 
cosmological parameters, it can also be considered as a probe of the clustering properties of the DM. Lensing dispersion has 



been investigated to search for the presence of compact objects in the Universe ( Linder, Schneider fc Wagoner 1988 



Ranch 
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1991; 



Metcalf . 



Silk 1999[ ^eljak fc Holz 1999 ). The possibility of determining the fraction of macroscopic DM using future 
samples of SNe la has also been explored (Vlortsell et al. 2001). SNAP should intensively observe SNe up to z ~ 1.7. In one 
year of study, this space-born mission should be able to discover ~ 2350 SNe, most of which in the region 0.5 ~ z ~ 1.2. The 
discrimination of models of Universe with different fractions of compact objects is mainly based on the shift in the peak of 
the lensing dispersion ( ^eljak fc Holz 1999 ; Mortsell et al. 2001): a shift of ~ 0.01 mag in the peak of the lensing dispersion in 
the projected SNAP sample towards lower amplifications corresponds to a growth of 20% in the fraction of macroscopic DM 
in a flat Universe with Qm = 0-3 and a cosmological constant (see figure (4) in Mortsell et al. (2001)). In Fig. we plot 
the dispersion in amplification, for the projected redshift distribution of SNe according to the SNAP proposal, in a Universe 
with Qm = 0.3 filled in with macroscopic DM . High de-amplification are preferred in the case of a cosmological constant, 
when the maximum of the distribution is depleted and the mode is shifted away from the mean with respect to dark energy 
with wq > —1. Changing from wq — —1 to wq — —1/2, the peak of the distribution moves for ~ 0.015 mag towards higher 
amplifications. So, a significant reduction in the fraction of compact object can be mimed by quintessence with wq > — 1. 
Since quintessence reduces the dispersion of gravitational lensing, it also reduces the ability to distinguish between microscopic 
and macroscopic DM from the shape of the amplification dispersion. Both quintessence and microscopic DM reduce the bias 
towards the empty beam value and the high magnification tail and their effect is of the same order. A Universe with an high 
fraction of macroscopic objects can be misleadingiy interpreted as one with dark energy with large negative pressure. 
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